Critical speeding-up in a local dynamics for the random-cluster model 
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We study the dynamic critical behavior of the local bond-update (Sweeny) dynamics for the 
Fortuin-Kasteleyn random-cluster model in dimensions d = 2, 3, by Monte Carlo simulation. We 
show that, for a suitable range of q values, the global observable S2 exhibits "critical speeding-up": 
it decorrelates well on time scales much less than one sweep, so that the integrated autocorrelation 
time tends to zero as the critical point is approached. We also show that the dynamic critical 
exponent z cxp is very close (possibly equal) to the rigorous lower bound a/u, and quite possibly 
smaller than the corresponding exponent for the Chayes-Machta-Swendsen-Wang cluster dynamics. 
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Dynamic processes in statistical mechanics typically 
undergo critical slowing- down the autocorrelation 
(relaxation) time r diverges as the critical point is ap- 
proached, most often like r ~ £ z , where £ is the spatial 
correlation length and z is a dynamic critical exponent. 
In this Letter we would like to draw attention to the con- 
verse (and quite unexpected, at least to us) phenomenon 
of critical speeding-up: some observables O can exhibit 
strong decorrelation on time scales much less than one 
sweep, so that the dominant relaxation modes equilibrate 
faster (in natural time units) near criticality. As a con- 
sequence, the integrated autocorrelation time Ti ntj o can 
in some cases tend to zero, so that the dynamic critical 
exponent z- lrA ,o is negative. These behaviors also have 
practical implications for the efficiency of Monte Carlo 
simulations [2|, |3[ near the critical point. 

More precisely, we shall exhibit these phenomena for 
the Fortuin-Kasteleyn (FK) random-cluster model 0, @ 
with a local bond- update dynamics The random- 
cluster model is a correlated bond-percolation model that 
is closely related to the Potts spin model @, S]- As such, 
it plays a major role in the theory of critical phenomena, 
especially in two dimensions where it arises in recent de- 
velopments of conformal field theory Q via its connection 
with stochastic Loewner evolution (SLE) [l(| 11 1 . 

The random-cluster model with parameters q, v > is 
defined on any finite graph G = (V, E) by the partition 
function 



Z = 



E 

ACE 



q KA) v \A\ 



(1) 



where A is the set of "occupied bonds" , \A\ is the number 
of occupied bonds, and k(A) is the number of connected 
components ("clusters") in the graph (V,A). For q — 1 



this reduces to independent bond percolation 12| with 
occupation probability p — v/(l + v): for integer q > 2 it 
provides a graphical representation [J] of the g-state fer- 
romagnetic Potts model with nearest-neighbor coupling 



positive values of q, integer or noninteger, to be studied 
within a unified framework. 

The simplest dynamics for the random-cluster model 
is the local bond-update dynamics first used by Sweeny 
[6[: choose a bond e S E at random, erase its current 
occupation state, and then give it a new occupation state 
according to the conditional distribution of (Q} with the 
other bonds held fixed [13(. In detail, this means that 
e will become occupied with probability v/(X + v) [resp. 
v/(q + v)] in case the endpoints of e are (resp. are not) 
already connected by a path of occupied bonds not using 
e. The efficient implementation of this connectivity check 
leads to nontrivial algorithmic questions [14j that we will 
discuss in detail elsewhere (l5| . 

In two dimensions, the behavior of the ferromagnetic 
Potts/random-cluster model is fairly well understood, 
thanks to a combination of exact solutions [l|l] , Coulomb- 
gas methods 17 1 and conformal field theory But 
in dimensions d > 3, many important aspects remain 
unclear, including the location of the crossover between 
second-order and first-order behavior [18| : the nature of 
the critical exponents and their dependence on q; the 
value of the upper critical dimension for noninteger q; and 
the qualitative behavior of the critical curve v c {q) near 
q = [l|| . Monte Carlo simulations using the Sweeny Q 
and Chayes-Machta [13] algorithms will likely play an 
important role in elucidating these problems. 

In this Letter we present the results of Monte Carlo 
simulations using the Sweeny dynamics, on d-dimensional 
simple hypercubic lattices of linear size L with periodic 
boundary conditions. We shall measure time in units of 
"hits" of a single bond, but we stress that the natural unit 
of time is one "sweep" of the lattice, consisting of dL d 
hits. For any observable O, we define the unnormalized 
autocorrelation function at time lag t, 



C (t) = (O s O s+t ) - {Of 



(2) 



J, where v 



O 0J 



1. The random-cluster model thus 



provides an extension of the Potts model that allows all 



where expectations are taken in the stationary stochastic 
process (i.e., in equilibrium), and the normalized autocor- 
relation function poit) — Co(t)/Co(0). We then define 



2 



the exponential autocorrelation time 



1*1 



and the integrated autocorrelation time 



Tint.O 



t— — oo 



Po{t) ■ 



(3) 



(4) 



Typically all observables O (except those that, for sym- 
metry reasons, are "orthogonal" to the slowest mode) 
have the same value T cxp .e> = r cxp . However, they may 
have very different amplitudes of "overlap" with this 
slowest mode; in particular, they may have very differ- 
ent values of the integrated autocorrelation time, which 
controls the efficiency of Monte Carlo simulations 
We define dynamic critical exponents z oxp and z- mtt o by 
T cx P ~ £ Zcxp and T intj o ~ £ Zi "t. , where time is measured 
in "sweeps" . On a finite lattice at criticality, £ can here 
be replaced by L. 

In our simulations we measured a variety of observ- 
ables, among which are the number of occupied bonds 
J\f = \A\ and the sum of squares of cluster sizes 52 = 
J2 \C\ 2 . It is well known 0] that (5 2 ) = V\, where 
V = L d is the volume and \ is the Potts-model suscepti- 
bility. A simple variational argument [2~0 ] shows that, in 
the Sweeny dynamics, T cxp > TintjAf > const x Ch, where 
Ch is the specific heat and time is measured in "sweeps" ; 
hence z cxp > z int ^ > a/u. 

We began by performing simulations on the square lat- 
tice (d — 2) at the exact critical point v c (q) = ^fq [l6| 
for q = 0.0005, 0.005, 0.05, 0.2, 0.5, 1, 1.5, 2, 2.5, 3, 3.5 and 
a variety of lattice sizes 4 < L < 1024 [22| . In all cases the 
autocorrelation function of Af is very close to a pure ex- 
ponential (Fig.[T]). The integrated autocorrelation times 
TintjA/" are shown as a function of q and L in Fig. [21 and the 
corresponding dynamic critical exponents z cxp « 2int,JV 
are shown in Table Q] The estimated exponents are only 
slightly larger than the lower bound a/v, and could con- 
ceivably be equal to it 24]. Perhaps surprisingly, these 
exponents are slightly smaller than those found recently 
[ia ] for the Chayes-Machta-Swendsen-Wang cluster 
algorithm. 

A more interesting and unusual dynamic behavior is 
exhibited by the observable £2- In Fig. [3] we plot the au- 
tocorrelation function ps 2 (t) for q = 0.2. Clearly ps 2 {t) 
exhibits a fast decay in a time much less than a single 
sweep (i.e. of order L w for some w < 2) as a prelude to 
the ultimate exponential decay e _t / Tcxp . To analyze this 
short-time behavior, we plot ps 2 (t) versus t/L w and ad- 
just the exponent w until all the points fall on a scaling 
curve ps 2 (t) — f{t/L w ) in the limit L — > oo. We find 
w ~ 0.99. Furthermore, the function / is very close to 
f(x) = (1 + ax)~ r with a = 0.55 and r = 1.42 (Fig. [J); 
in any case it behaves like f(x) ~ x~ r as x — > oo. 
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Figure 1: Autocorrelation function pjv(t) versus t/L 2 for 
the critical two-dimensional random-cluster model at q = 0.2, 
where time t is measured in "hits" . 
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Figure 2: Integrated autocorrelation times Ti n t,AT for the crit- 
ical two-dimensional random-cluster model, as a function of q 



and L. For q = 1 the analytical result is n n t,Af 
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0.0005 





-1.9576 


0.77 4.83 


-1.23 


1.2376 


1.9965 
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0.79 


4.18 


-1.21 


1.2111 


1.9891 
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2.84 


-1.12 
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1.42 


-1.01 


1.0168 


1.9417 
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-0.8778 


1.11 
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-0.71 


0.8904 


1.9172 


1.0 





-0.5000 


1.26 


0.43 


-0.32 


0.7500 


1.8958 


1.5 





-0.2266 


1.36 


0.25 


-0.16 


0.6398 


1.8832 
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(log) 


(log) 


1.49 


0.15 


-0.08 


0.5417 


1.8750 


2.5 


0.26(1) 


0.2036 


1.64 


0.10 


0.20 


0.4474 


1.8697 


3.0 


0.45(1) 


0.4000 


1.84 


0.06 


0.41 


0.3500 


1.8667 


3.5 


0.636(2) 


0.6101 


2.04 


0.04 


0.61 


0.2375 


1.8662 



Table I: Estimated dynamic critical exponents for the two- 
dimensional random-cluster model as a function of q. Specific- 
heat exponent a/v and red-bond (resp. whole-cluster) fractal 
dimension d r cd (resp. dp) are shown for comparison [23| . 




Figure 3: Autocorrelation function ps 2 (t) versus t/L 2 for 
the critical two-dimensional random-cluster model at q = 0.2. 
The autocorrelation function p^r(t) is shown for comparison. 
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Figure 4: Autocorrelation function ps 2 (t) versus l+at/L w for 
the critical two-dimensional random-cluster model at q = 0.2, 
with w = 0.99 and a = 0.55. 



Finally we can analyze the universal crossover from 
short-time to long-time behavior, which we hypothesize 
is of the form p Sz (t) = f(t/L w )g(t/L d+z ^) 1 by plotting 
(l+at/L w ) r p S2 (t) versus t/L d+Zcx p. A fairly clear scaling 
curve is seen (Fig. [5]), though it is noisy for large lattices. 
Using this scaling Ansatz to compute the area under the 
curve of ps 2 (t) , we conclude that 



^int,<S2 



r(w - d) + (1 - r)z cxp if r < 1 
w — d if r > 1 



(5) 



Similar analyses for the other values of q yield the expo- 
nents reported in Table fl] Note that critical speeding-up 
(w < d) and critical slowing-down (z oxp > 0) can coexist. 

Critical speeding-up also occurs in the random-cluster 
model in dimensions d > 2. We simulated the three- 



dimensional random-cluster model 22| at the estimated 
critical points w = v/q = 0.43365 for q = [l9j . p 



v/(l +v) = 0.2488126 for q 
with P = 0.22165455 for q = 




and v 



1 



t/L 

Figure 5: (1 + t/L w ) r ps 2 (t) versus t/L 2 for the critical two- 
dimensional random-cluster model at q = 0.2, with w — 0.99, 
a — 0.55 and r = 1.42. The autocorrelation function pj^(t) is 
shown for comparison. 
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1.52 1.04 
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0.35(1) 
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2.55 0.08 


0.29 


0.757(2) 2.4816(1) 



using lattice sizes 



Table II: Same as Table |TJ for the three-dimensional random- 
cluster model, ah/ and d F from ^ [UHil for q = 0, 1, 2; 
drcd from [2|| [20] for q = 1,2. 



4 < L < 256 for g = 0, 1 and 4 < L < 64 for q = 2. The 
qualitative behavior was the same as in d = 2, and the 
estimated exponents are shown in Table [ill Our value 
of z CKp for q — 2 is consistent with that of Wang, Kozan 
and Swendsen 0]. In fact, T mt ^/CH for q = 2 is close 
to constant, so it is conceivable that z cxp = a/V exactly. 

In retrospect it is not surprising that a "global" ob- 
servable like 1S2 could exhibit significant decorrelation in 
a time much less than one sweep. After all, FK clus- 
ters are fractals: a large cluster can sometimes be broken 
into two large pieces by one or a few bond deletions, and 
two large clusters can sometimes be joined by one or a 
few bond insertions. This reasoning suggests correctly 
that the critical speeding-up should be strongest when 
the cluster is most fragile, i.e. the red-bond fractal di- 
mension d Te d [23| is largest, namely at small q. 

We can pursue this idea further and suggest that the 
decorrelation of 1S2 is caused principally by hitting a few 
(order 1) red bonds: this takes a time oc L d ~ dicd , so we 
predict w — d — c? rc d- Our data (Tables HI and Hlf are in 
amazingly good agreement with this prediction for q < 2, 
i.e. when z cxp = 0. However, they deviate from it when 
Zcxp > 0, for reasons that we do not yet understand. Note 
in particular that for d > 6 and q = 1 (resp. < q < 1), 
we expect <i re d = Vto — 2 (resp. d Ic a > 2) and hence 
w = d — 2 (resp. w < d — 2). 

We lack, at present, any theory (or even any numerol- 
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